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^h , Abstract. Analytical arguments and dynamic Monte Carlo simulations show that 

the microscopic structure of field-driven Solid-on-Solid interfaces depends strongly on 

the details of the dynamics. For nonconservative dynamics with transition rates that 

factorize into parts dependent only on the changes in interaction energy and field 

energy, respectively (soft dynamics), the intrinsic interface width is field-independent. 

For non-factorizing rates, such as the standard Glauber and Metropolis algorithms 

(hard dynamics), it increases with the field. Consequences for the interface velocity 

and its anisotropy are discussed. 
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Surfaces and interfaces moving under far-from-equilibrium conditions are important in 
the formation of patterns and structures. For example, domain growth by the motion of 
defects such as grain boundaries or dislocations influences the mechanical properties 
of metals JI| and the structures formed during phase transformations in adsorbate 
systems |2| , while the propagation of domain walls influences the switching dynamics of 
magnetic nanoparticles and ultrathin films ||. Recently, interface dynamics was even 
used to analyze the scalability of discrete-event simulations on parallel computers @]. 
The importance of moving interfaces in a wide variety of fields has inspired enormous 
interest in their structure and dynamics ||. However, despite the fact that many 
interface properties, such as mobility and chemical activity, are determined by the 
microscopic structure, interest has mostly focused on the large-scale interface structure 
and its universal properties. In this Letter we show that the microstructure of a 
moving interface can be dramatically influenced by seemingly minor modifications of 
the growth mechanism, with measurable consequences for such macroscopic quantities 
as the interface velocity and its anisotropy. 

The detailed microscopic mechanism of the interface motion is usually not known, 
and it is therefore common to mimic the essential features in a dynamic Monte Carlo 
(MC) simulation of a model stochastic process |]]. In doing so, there are two important 
distinctions that must be made. One is between those transition probabilities that only 
depend on the energy difference between the initial and final states (often referred to as 
Metropolis, Glauber, or heat-bath rates || 0, |J), and those that involve an activation 
barrier between the two states (often referred to as Arrhenius rates) 0, |j, |8|, |[J . Another 
distinction is between dynamics that do not conserve the order parameter, such as the 
Metropolis and Glauber single-spin flip algorithms, and conservative dynamics, such as 
the Kawasaki spin-exchange dynamic Rj . Once it is decided to which of these categories 
the system belongs, the dynamic is often chosen on the basis of convenience. 

It is well known that different microscopic dynamics can yield different equilibration 
paths and equilibrium fluctuations [|TUJ (cluster vs local MC algorithms being the most 
extreme example 0) and even noticeable differences in the steady-state microstructure 
II ||. Nevertheless, the general expectation is that, if no additional parameters 
(such as an activation barrier or a diffusion rate) is introduced into the physical 
model, observables are only affected quantitatively. However, here we demonstrate 
striking qualitative differences between the microstructural consequences of two different 
stochastic dynamics that do not involve a transition barrier, even though both obey 
detailed balance and neither involves order-parameter conserving moves. For this 



demonstration we use an unrestricted Solid-on-Solid (SOS) interface |TT[ separating 
uniform spin-up and spin-down phases in a ferromagnetic S = 1/2 Ising system 
on a square lattice of unit lattice constant. The interface consists of integer- valued 
steps parallel to the y-direction, 5(x) e {— oo, +oo}, where the step position x is 
also an integer. Its energy is given by the nearest-neighbour Ising Hamiltonian with 
anisotropic ferromagnetic interactions, J x and J y in the x and y direction, respectively: 
H = ~J2x,y s x,y(JxS x +i,y + JyS x ,y+i + H). The two states at site (x, y) are denoted 
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by the spin s x>y = ±1, J2x,y runs over a ^ sites, and H is the applied field. We 
take s x>y = —1 on the side of the interface corresponding to large positive y, so 
that the interface on average moves in the positive y direction for H > 0. In the 
equivalent lattice-gas language, s x>y = +1 corresponds to solid and s x>y = — 1 to gas or 
liquid, and H is proportional to a chemical-potential imbalance. The interface evolves 
through single spin flips (in lattic-gas language: adsorption/desorption events) that 
occur with probability Wj/^A-E 1 ] where AE is the energy change that would result 
from the transition, constrained by the requirement that the probabilities obey detailed 



balance (see fT2|, [13[ for details). To restrict the accessible configurations to a simple 
SOS interface without bubbles or overhangs, transitions are allowed only for spins that 
have one single broken bond in the ^/-direction. Regardless of the details of the transition 
probabilities, which lead to dramatic microscopic structural differences as demonstrated 
below, the macroscopic structure of the interfaces belongs to the Kardar-Parisi-Zhang 



KPZ) dynamic universality class |5|, 14 1. We emphasize that the absence of energy 



barriers and order-parameter conserving diffusion (spin-exchange) moves sets our model 
clearly apart from commonly studied models of Molecular Beam Epitaxy (MBE) |], ||. 
It is closely similar to the diffusion-free dynamical model studied in three dimensions in 
Rcf. [|rj| which, however, does not consider the microscopic interface structure in detail. 



Recently |T2| we introduced a mean-field approximation for the driven-interface 
microstructure, in which the probability density function (pdf ) for the height of a single 



step takes the Burton-Cabrera- Frank (BCF) form [ |TTf , 

p[S(x)} = Z- l X^ x) \ e^ mx) , (1) 

where the Lagrange multiplier 7(0) enforces (S(x)) = tan0 independent of x, 
corresponding to an overall angle between the interface and the x-axis. The width 
parameter X, which for a driven interface can depend on both H and the temperature 
T, is discussed below. The partition function in ([!]) is 

+00 1 y2 

Z - V X |5| p 7W<5 - - - (2) 

~ ^ A 6 ~ 1-2X cosh 7 (0)+ X 2 [) 



6 



with 



7(0) __ (1 + X 2 ) tan 0+V(l-X 2 ) 2 tan 2 + 4X 2 
6 2X(l + tan0) ' [ } 

For = 0, 7(0) = 0, yielding Z(0) = (1 + X)/(l - X). 

In the approximation of jl^l , individual steps are assumed to be statistically 
independent (like in the original BCF model for equilibrium interfaces). The width 
parameter X is then obtained self-consistently as [113, [L3| 



X(T m Y m [ e-^W[(3(-2H - U x )] + e^WiP(+2H - U x )} \ 1/2 
X(T,tf)-X (T)j Wm -2H-4J X )] + Wm+2H-U x )] j ' (4) 

where X (T) = e~ 2/3Jx with (3 = l/k^T (k-Q is Boltzmann's constant) is the BCF zero- 
field equilibrium value of X. For H = 0, X(T, H) reduces to Xq(T). For instance, using 
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the Glauber transition probability, 



e - P AE 



we get 



W G (S X ,y -> S X ,y) = - ~ e _ /3A£ , (5) 

Y f T ff ^ y m / e 2/3J - cosh(2/?ff) + e~^ 1 1/2 

Equation (0) exemplifies a class of dynamics known as /iani in the literature on 
nonequilibrium lattice models [IE ], in which the transition probabilities depend directly 



on the total energy change, AE. In a different class, known as soft, the probabilities 
factorize into a part due only to the change in the field energy, AEh oc H, and a part 
due only to the change in the interaction energy, AEj oc J x . One example is the soft 
Glauber dynamic, 

e -(3AE H e -PAEj 

W SG (s x , y -> -s x ,„) = - - e _ pAEH ■ x - e _ pAEj , (7) 

which also obeys detailed balance. 

The important point is that when a factorizing transition probability such as (|7|) 
is used in ([|), all dependence on H cancels due to the detailed balance. Whether one 
uses a soft or a hard dynamic thus strongly affects the intrinsic interface width and 
properties that depend on it, such as the propagation velocity. The need to use a soft 
dynamic in cases where the "field" represents a chemical-potential difference has been 
recognized in some MC studies of crystal growth |T7], [18], [L9| . 



Next we illustrate the large differences between hard and soft dynamics by 
comparing analytic approximations and dynamic MC simulations of driven SOS 
interfaces with J x =J y =J, using the hard and soft Glauber dynamics, respectively. In 
particular, we consider the field- and temperature dependencies of the intrinsic interface 
width, represented by the average absolute value of the step height, and the resulting 
interface velocity and its anisotropy. 

In figure [I] we show numerical evidence supporting our predictions for the average 
step height, which for (ft = is related to X as (\5\) = 2X/(1 — X 2 ), under both the 
hard and soft Glauber dynamics at two different temperatures. (Simulational details 
are given below.) The difference between the two dynamics is striking: the step heights 
for the soft dynamic are independent of H, in contrast to the strong iJ-dependence 
produced by the hard dynamic. Additional confirmation of the functional form of the 
single-step pdf, ([[]), is obtained from the MC data by calculating (\8\), both directly 
by summation over the numerically obtained pdf, and from the probability of zero step 
height for = as (\S\) = {p[0] _1 — p[0]}/2. This result is obtained by combining 
p[0] = Z^ 1 with the above relation for (\S\) in terms of X. 

The spins in an Ising system can be classified by the value of AE, which is uniquely 
determined by s and the number of broken x- and y-bonds. The active spin classes all 
have one broken y-bond and can be labeled by s=±l and the number of broken x- 
bonds, j G {0, 1, 2}, as js. The classes and corresponding AE are given in the first two 
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columns of table |1|. The details of our "active-site" implementation of the spin-class 
based discrete-time ra-fold way algorithm |K]] are given in |I2| , the only difference here 
being the restriction to flips only at sites with one broken ?/-bond. The data shown here 
are for interfaces of size L x =10 4 in the x-direction, which were allowed 5 000 n-fold way 
updates per active site to reach stationarity, after which pdfs, class populations, and 
velocities were collected over 50 000 updates per active site. 

Assuming up-down symmetry of the interface, the mean perpendicular propagation 
velocity is given by 

(v ± (T,H,4>)} = co S( j>J2(n(js)}(v y {j)} (8) 

3 

in units of inverse MC steps per site (MCSS -1 ). Here, the average spin-class populations, 
(n(js)), are calculated from the single-step pdf, assuming statistical independence 
|12|] , and the results are given in the third and fourth columns of table |l|. The 



corresponding contributions to the interface velocity in the y-direction, (v y (j)) = 

{W [(3AE(j-)} —W\ftAE(j+)]}, are given in the fifth and sixth column of table |l] 
for the hard and soft Glauber dynamic, respectively. The special case of 4> = leads to 
compact formulas: 

'^°»'''wj" + i t — -" ' (9) 

with X from @ for the hard Glauber dynamic, and 

MT,#,0)) SG = tanhOStf)^^ ( 10 ) 



smh(2f3J x ) 
COBh(/3H) 



1 + X n 2 



for the soft Glauber dynamic. These results are compared with the corresponding MC 
data in figure ^|. The difference between the two dynamics is clear. 

The angular dependencies of (v±) for the two dynamics are illustrated in figure §. 
The theoretical results represent (|j) with terms from table [3]. The hard dynamic yields 
an anisotropy which changes from that of the single-step model for weak H, to the 



reverse anisotropy of the Eden model |H], |22[ for strong H. There is no such change 



for the soft dynamic, which shows a single-step like anisotropy for all fields. Effects 
of different transition probabilities on nonequilibrium growth shapes have already been 
observed in systems with conserved order parameter |T(J. It is likely that the different 
anisotropies resulting from the hard and soft dynamics will lead to different growth 
shapes in the case of nonconserved order parameter, as well. 

With hard dynamics, the up-down symmetry of the interface is gradually destroyed 



with increasing field j|, [T2], [13], [£J, leading to different class populations for s=+l and 
s=— 1. Our MC simulations show this skewness to be absent for the soft dynamic. 
This makes us believe that the field dependence in (^) is exact in the following sense. 
If the interface has up-down symmetry and is driven with a soft dynamic, then the 
generated interface shapes are independent of H. Consequently, the only dependence of 
the velocity on H is in the difference between the spin-flip probabilities of a down spin 
and an up spin, which yields ta.nh(/3H ) for the Glauber dynamic. 
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In summary, we have presented analytical and dynamic MC results that show 
remarkable qualitative differences in the microscopic structure of field-driven SOS 
interfaces in the KPZ dynamic universality class, depending on whether the local 
stochastic spin dynamic is hard (such as the standard Glauber and Metropolis dynamics) 
or soft. In the former case, (\S\) depends strongly on H, while in the latter it has no in- 
dependence at all. We emphasize that the difference between the hard and soft dynamic 
is purely dynamical and does not involve introduction into the physical model of any 
new parameters, such as an activation barrier or a diffusion rate. If microstructural 
information is desired from a dynamic MC simulation, great care must therefore be 
exercised, both in choosing the transition probabilities and in interpreting the results. 



Analogous effects for other models, such as Ising interfaces [12], and in three dimensions 



and the possibility that the resulting growth shapes may be different [HJ, are left for 



future study. Further details on the microstructure of the SOS interface under the hard 



Glauber dynamic will be reported elsewhere |L3 . 



Acknowledgments 

We thank G Brown, G Buendia, S J Mitchell, M A Novotny, and K Park for useful 
comments. PAR appreciates hospitality at Kyoto University. Supported in part by 
US NSF Grant No. DMR-9981815, and by Florida State University through MARTECH 
and CSIT. 

References 



[1] 
[2] 
[3 



[4] 
[5] 

m 

[7] 

m 

[9] 
[10 

[11 

[12 
[13 
[14 
[15 



Mendelev M I and Srolovitz D J 2000 Acta Mater. 48 3711; 2001 ibid 49 589; 2001 ibid 49 2843 
Mitchell S J, Brown G and Rikvold P A 2001 Surf. Set. 471 125 

Rikvold P A, Brown G, Mitchell S J and Novotny M A in press Nano structured Megnetic Materials 
and their Applications ed D Shi, B Aktas, L Pust and F Mikailov (Berlin: Springer) (preprint 



cond-mat/0110103D 



Korniss G, Toroczkai Z, Novotny M A and Rikvold P A 2000 Phys. Rev. Lett. 84 1351 
Barabasi A-L and Stanley H E 1995 Fractal Concepts in Surface Growth (Cambridge: Cambridge 

Univ. Press) and references therein 
Siegert M and Plischke M 1994 Phys. Rev. E 50 917 
Landau D P and Binder K 2000 Monte Carlo Simulations in Statistical Physics (Cambridge: 

Cambridge Univ. Press) 
Shim Y and Landau D P 2001 Phys. Rev. E 64 036110 
Kang H C and Weinberg W H 1989 J. Chem. Phys. 90 2824; Fichthorn K A and Weinberg W H 

ibid 95 1090 
Shochet O, Kassner K, Ben-Jacob E, Lipson S B and Muller-Krumbhaar H 1992 Physica A 181 

136 and references therein 
Burton W K, Cabrera N and Frank F C 1951 Phil. Trans. Roy. Soc. (London) Ser. A 243 299 
Rikvold P A and Kolesik M 2000 J. Stat. Phys. 100 377 
Rikvold P A and Kolesik M in preparation 

Kardar M, Parisi G and Zhang Y-Z 1986 Phys. Rev. Lett. 56 889 
Jiang Z and Ebner C 1992 Phys. Rev. B 45 6163 



Letter to the Editor 7 

[16] Marro J and Dickman R 1999 Nonequilibrium Phase Transitions in Lattice Models (Cambridge: 

Cambridge Univ. Press) ch 7 
[17] Clio H, Grossmann B and Grant M 1990 Phys. Rev. Lett. 64 1262 
[18] Kotrla M and Levi A C 1991 J. Stat. Phys. 64 579 
[19] Hontinfinde F, Krug J and Touzani M 1997 Physica A 237 363 
[20] Bortz A B, Kalos M H and Lebowitz J L 1975 J. Comput. Phys. 17 10 
[21] Meakin P, Jullien R and Botet R 1986 Europhys. Lett. 1 609 
[22] Hirsch R and Wolf D E 1986 J. Phys. A 19 L251 
[23] Neergaard J and den Nijs M 1997 J. Phys. A 30 1935 



Table 1. The three spin classes that contribute to the interface velocity in the SOS 
model (column 1), together with the corresponding energy changes resulting from a 
successful spin flip (column 2, with upper sign corresponding to initial s = — 1 and 
lower sign to s = +1), the average class populations for general tilt angle 4> (column 
3) and for <fi ~ (column 4), and the contributions to the interface velocity in the y- 
direction for the hard Glauber dynamic (column 5) and for the soft Glauber dynamic 
(column 6). In columns 3 and 4, X corresponds to the width parameter X(T,H) for 
the specific dynamic used, given by (0) in the general case, (0) for the hard Glauber 
dynamic, and Aq(T) = e~ 2)3Jx for soft dynamics. 
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Figure 1. Average step height (\5\) vs H for 0=0 at T=0.2T C and 0.6T C , where T c 
is the exact Ising critical temperature. Solid curves, obtained as (\6\) = 2X/(1 — X 2 ) 
with X — Xq(T,H), are theoretical predictions for the hard dynamic. Dashed lines, 
similarly obtained with X = Xq(T), are theoretical predictions for the soft dynamic. 
The MC data are from direct summation over the simulated single-step pdfs (filled 
symbols) and from p[0] using (\5\) = {p[0] _1 — p[0]} /2 (empty symbols). In this and 
the following figures, the statistical uncertainty is much smaller than the symbol size. 




Figure 2. The perpendicular interface velocity (v±) vs H for <fi = at T=Q.2T C 
and 0.6T C . Theoretical results for the hard dynamic are from (||) (solid curves), with 
MC data shown as squares. Theoretical results for the soft dynamic are from (|hJ) 
(dashed curves), with MC data shown as diamonds. 
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Figure 3. Angle dependence of the perpendicular interface velocity (v±) at T=0.2T C 
for (from below to above) H/J = 0.1, 0.5, 1.0, 1.5, 2.0, and 2.5. Curves are theoretical 
predictions, data points MC results, (a) Hard dynamic, (b) Soft dynamic. On this 
scale, the soft-dynamic results overlap for H/J > 1.5. 



